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Abstract 


We discuss the reconstruction of piecewise smooth data from its (pseudo-) spectral informa¬ 
tion. Spectral projections enjoy superior resolution provided the data is globally smooth, while 
the presence of jump discontinuities is responsible for spurious 0(1) Gibbs oscillations in the 
neighborhood of edges and an overall deterioration to the unacceptable first-order convergence 
rate. The purpose is to regain the superior accuracy in the piecewise smooth case, and this is 
achieved by mollification. 

Here we utilize a modified version of the two-parameter family of spectral mollifiers intro¬ 


duced by Gottlieb & Tadmor |GoTa85|. The ubiquitous one-parameter, finite-order mollifiers 
are based on dilation. In contrast, our mollifiers achieve their high resolution by an intricate 
process of high-order cancelation. To this end, we first implement a localization step using edge 
detection procedure, |GeTa00a, GeTaOOb . The accurate recovery of piecewise smooth data is 


then carried out in the direction of smoothness away from the edges, and adaptivity is responsible 
for the high resolution. The resulting adaptive modifier greatly accelerates the convergence rate, 
recovering piecewise an alytic da ta within exponential accuracy while removing spurious oscilla¬ 
tions that remained in [ GoTa85|| . Thus, these adaptive mollifiers offer a robust, general-purpose 
“black box” procedure for accurate post processing of piecewise smooth data. 
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1 Introduction 

We study a new procedure for high resolution recovery of piecewise smooth data from its (pseudo- 
spectral information. The purpose is to overcome the low-order accuracy and spurious oscillations 
associated with Gibbs phenomena, and to regain the superior accuracy encoded in the global 
spectral coefficients. 

A standard approach for removing spurious oscillations is based on mollification over a local region 
of smoothness. To this end one employs a one-parameter family of dilated unit mass modifiers of 
form (fig = tp(x/9)/9. In general, such compactly supported modifiers are restricted to finite-order 
accuracy, \(pg * f(x) — f(x )| < C r 6 r , depending on the number r of vanishing moments tp has. 
Convergence is guaranteed by letting the dilation parameter, 9 { 0. 

In |PoTa85| we introduced a two-parameter family of spectral modifiers of the form 


VwO®) = \p(^) d p(p)- 


9' 


Here p(-) is an arbitrary Cq°(— 7r, 7r) function which localizes the p-degree Dirichlet kernel D p (y ) := 
• The parameter — the dilation parameter 9 need not be small in this case, in fact 
9 = 9{x) is made as large as possible while maintaining the smoothness of p(x — 9■)/{■). Instead, it 
is the second parameter - the degree p, which allows the high accuracy recovery of piecewise smooth 
data from its (pseudo-)spectral projection, P^f{x). The high accuracy recovery is achieved here 
by choosing large p’s, enforcing an intricate process of cancelation as an alternative to the usual 
finite-order accurate process of localization. 

In §| we begin by revisiting the convergence analysis of | GoTa.85|1 . Spectral accuracy is achieved 
by choosing an increasing p y/N, so that 'ipp.o has essentially vanishing moments ad orders, 
f y s i’ P ,e(y)dy = 5 s o + C s ■ N s / 2 , Vs, yielding the ’infinite-order’ accuracy bound in the sense of 
\^ P ,e* PnHx) - f(x)\ <C S -N~ S / 2 , Vs. 


Although the last estimate yields the desired spectral convergence rate sought for in |GoTa85| l , it 
suffers as an over-pessimistic restriction since its derivation ignores the possible dependence of p on 
the degree of local smoothness, s, and the support of local smoothness, ~ 9 = 9(x). In §|| we begin 
a detailed study on the optimal choice of the (p, 9) parameters of the spectral modifiers, Vv,0 : 

• Letting d(x) denote the distance to the nearest edge, we first set 9 = 9{x) ~ d(x) so that 
-k Pn f(x) incorporates the largest smooth neighborhood around x. To find the distance 
to the nearest discontinuity we utilize a general edge detection procedure, |GeTa9S|, peTaOOa, 


GeTaOObj , where the location (and amplitudes) of ad edges are found in one global sweep. 


Once the edges are located it is a straightforward matter to evaluate, at every x, the appro¬ 
priate spectral parameter, 9(x) = d(x)/ir. 

Next, we turn to examine the degree p, which is responsible for the overall high accuracy by 
enforcing an intricate cancelation. A careful analysis carried out in §L1 leads to an optimal 
choice of an adaptive degree of order p = p{x) ~ d(x)N. Indeed, numerical experiments 


reported back in the original |GoTa85|] and additional tests carried out in §L2 below and 


which motivate the present study, clearly indicate a superior convergence up to the immediate 
vicinity of the interior edges with an adaptive degree of the optimal order p = p{x) ~ d{x)N. 


Given the spectral projection of a piecewise analytic function, Sj\rf(-), our 2-parameter family of 
adaptive modifiers, equipped with the optimal parameterization outlined above yields - consult 


Theorem 3.1 below 
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H’pfi * S^f^x) — f(x) | < Const ■ d{x)N ■ e pV d ( x l N . 


The last error bound shows that the adaptive modifier is exponentially accurate at all x’s except for 
the immediate 0(l/A)-neighborhood of the jumps of /(•) where d(x) ~ l/N. We note in passing 
the rather remarkable dependence of this error estimate on the Cfif regularity of p(-). Specifically, 
the exponential convergence rate of a fractional power is related to the Gevrey regularity of the 
localizer /?(•); in this paper we use the GVregular cut-off p c (x) = exp(cx 2 /(x 2 — 7r 2 )) which led to 
the fractional power 1/2. 

Similar results holds in the discrete case. Indeed, in this case, one can bypass the discrete Fourier co¬ 
efficients: expressed in terms of the given equidistant discrete values, {f(yv)}, of piecewise analytic 


/, we have - consult Theorem 3.2 below, 


7T 

'N 


2N-1 


^pA x - Vv)f{yv) - f( x ) I < Const ■ {d{x)N) 2 ■ e W d C) N _ 


v =o 


Thus, the discrete convolution 'Yh v if Pl o{x — y v ) f {y v ) forms an exponentially accurate near-by inter- 
polant0, which serves as an effective tool to reconstruct the intermediate values of piecewise smooth 
data. These near-by “expolants” are reminicient of quasi-interpolants, e.g., [BL92], with the em¬ 
phasize given here to nonlinear adaptive recovery which is based on global regions of smoothness. 


What happens in the immediate, 0(l/iV)-neighborhood of the jumps? in §[| we complete our 
study of the adaptive mollifiers by introducing a novel procedure of normalization. Here we enforce 
the first few moments of the spectral modifier, fjp = pD p to vanish, so that we regain polynomial 
accuracy in the immediate neighborhood of the jump. Taking advantage of the freedom in choosing 
the localizer /?(•), we show how to modify p to regain the local accuracy by enforcing finitely 
many vanishing moments of if p = pD p , while retaining the same overall exponential outside the 
immediate vicinity of the jumps. By appropriate normalization, the localized Dirichlet kernel we 
introduce maintains at least second order convergence up to the discontinuity. Increasingly higher 
orders of accuracy can be worked out as we move further away from these jumps and eventually 
turning into the exponentially accurate regime indicate earlier. In summary, the spectral modifier 
amounts to a variable order recovery procedure adapted to the number of cells from the jump 
discontinuities, which is reminiscent of the variable order, Essentially Non Oscillatory piecewise 
polynomial reconstruction in [ HEOC85 |. The currrent procedure is also reminicient of the h-p 
methods of Babuska and his collabrators, with the emphasize given here to an increasing number 
of global moments (p) without the u h” -refinement. The numerical experiments reported in §|3.2| and 


§ 4.3 confirm the superior high resolution of the spectral modifier if p p equipped with the proposed 
optimal parameterization. 


Acknowledgment. Research was supported in part by ONR Grant No. 
NSF grant #DMS01-07428. 
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2 Spectral Mollifiers 

2.1 The two-parameter spectral mollifier if Py o 

The Fourier projection of a 27r-periodic function /(•), 


1 Called expolant for short 
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SnHx) := y kX i A == ^ £ /(*)< 


,—ikx 


dx 


( 2 . 1 ) 


|fc|<JV 


enjoys the well known spectral convergence rate, that is, the convergence rate is as rapid as the 
global smoothness of /(•) permits in the sense that for any s we have/ 


\SnI(x) - f(x )| < Consi||/||c s 


1 


Vs. 


AT *- 1 

Equivalently, this can be expressed in terms of the usual Dirichlet Kernel 


Dn W := i E 


N 


2ir 


e ikx = sin(iv + l/2)a 
27rsin(x/2) 


( 2 . 2 ) 


(2.3) 


k=—N 

where S^f = Djy * /, and the spectral convergence statement in 

I D n * f(x) - f{x )| < Const||/1|c s • 

Furthermore, if /(•) is analytic with analyticity strip of width 2g, then Swf(x ) is characterized by 
an exponential convergence rate, e.g., 


recast into the form 


Vs. 


(2.4) 


Ta94] 


\Sn f{x) — f(x)\ < Const^ ■ Ne Nv . (2-5) 

If, on the other hand, /(•) experiences a simple jump discontinuity, say at xo, then Sjsrf(x) suffers 
from the well known Gibbs’ phenomena, where the uniform convergence of Swf(x) is lost in the 
neighborhood of xq, and moreover, the global convergence rate of SnI(x) deteriorates to first order. 
To accelerate the slow convergence rate, we focus our attention on the classical process of mollifi¬ 
cation. Standard modifiers are based on a one-parameter family of dilated unit mass functions of 
the form 


¥>«(*) : = ( 2 - 6 ) 

which induce convergence by letting 0 to zero. In general, \ipg * f(x) — f(x)\ < C r 6 r describes the 
convergence rate of finite order r, where (p possesses r vanishing moments 


J y s v{y)dy = b sQ s = 0,1,2,... , r - 1. 


(2.7) 


In the present context of recovering spectral convergence, however, we follow Gottlieb and Tadmor, 
|]GoTa85 |, using a two-parameter family of modifiers, fi P! g(x), where 6 is a dilation parameter, 
^pfi^x) = fi’pix / 9) / 6 , and p stipulates how closely fi P: g(x) possesses near vanishing moments. To 
form ip p (x), we let p(x) be an arbitrary C^° function supported in (— tt, n) and we consider the 
localized Dirichlet kernel 


fip(x) := p(x)D p (x). (2.8) 

Our two-parameter modifier is then given by the dilated family of such localized Dirichlet kernels 



2 Here and below we denote the usual ||/||c s := II/^IIl 00 - 
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According to (2.9), ijj p g consists of two ingredients, p{x) and D p (x), each has essentially separate 
role associated with the two independent parameters 9 and p. The role of p (|) is, through its 9- 
dependence, to localize the support of il> Pt g(x) to (~9tt,9tt). The Dirichlet kernel D p (x ) is charged, 
by varying p, with controlling the increasing number of near vanishing moments of i/> p g, and hence 
the overall superior accuracy of our mollifier. Indeed, by imposing the normalization of 


P(0) = 1, (2.10) 

we find that an increasing number of moments of il> P: g are of the vanishing order 0(p~^ s ~ 1 ' > ), 

/ it 6 rir 

y s ^p,e(y) d y= (yO) s p{y) D P (y)dy = D p -k(y9) s p(y)\ y=0 = 5 s0 + C s ■ p- {s ~ 1 '> Vs, (2.11) 

-irQ J — 7r 


where according to {2A), C s = C onst\\(y9) s p(y)\\c s ■ We shall get into a detailed convergence 
analysis in the discussion below. 

We conclude this section by highlighting the contrast between the standard, polynomially accurate 
mollifier, (^7j) and the spectral mollifiers (2.S). The former depends on one dilation parameter, 
9, which is charged of inducing a fixed order of accuracy by letting 9 | 0. Thus, in this case 
convergence is enforced by localization , which is inherently limited to a fixed polynomial order. 
The spectral mollifier, however, has the advantage of employing two free parameters: the dilation 
parameter 9 which need not be small - in fact, 9 is made as large as possible while maintaining 
p(x — 9y)f(y) free of discontinuities; the need for this desired smoothness will be made more evident 
in the next section. It is the second parameter, p, which is in charge of enforcing the high accuracy 
by letting p j" oo. Here, convergence is enforced by a delicate process of cancelation which will 
enable us to derive, in §0, exponential convergence. 


2.2 Error analysis for spectral mollifier 

We now turn to consider the error of our mollification procedure, E(N,p,9; f(x)), at an arbitrary 
fixed point x E [0, 27r) 


E{N,p,9;f(x))=E(N,p,9) := ^ p ,g * S N f(x) - f(x), (2.12) 

where we highlight the dependence on three free parameters at our disposal - the degree of the 
projection, N, the support of our mollifier, 9, and the degree with which we approximate an arbitrary 
number of vanishing moments, p. The dependence on the degree of piecewise smoothness of /(•) 
will play a secondary role in the choice of these parameters. 

We begin by decomposing the error into the three terms 

E(N,p, 9) = (f-k if pfi - /) + (S N f - f) -k (V'p.e - S N ^ P)d ) + ( S N f - f) * S N ^ p>e . (2.13) 

The last term, (Sjv/ — /) * -SVV’p.e, vanishes by orthogonality, and hence we are left with the first 
and second terms, which we refer to as the Regularization and Truncation errors, respectively 


E(N,p, 9) = (f-k -ip^g - /) + (S N f - f) -k (ij; Pt g - S N ^ P) g) 


R(N,p, 9) + T(N,p, 9). 


(2.14) 


Sharp error bounds for the regularization and truncation errors were originally derived in [ GoTa85 ] 
and a short re-derivation now follows. 

For the regularization error we consider the function 
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9x{y) ■= f{x - 9y)p{y) - f(x) 


(2.15) 


where f(x ) is the fixed point value to be recovered through mollification. Applying (2.4) to g x {- 
while noting that g x ( 0) = 0, then the regularization error does not exceed 


\R(N,p : 9)\ := |/*^ P) e-/| = 


[ [fix - 9y)piy ) - f(x)]D p (y)dy 

J —7V 


1 


= \D p -k g x iy)\ y=0 \ = \iS p g x iy ) - p. T (p))| y= o| < C'onst||^(y)||c= • — r - (2.16) 


Applying Leibnitz rule to g x iy ), 




laLtell <£( * ) o k \f fk) (x - 0v)\ ■ |p ( -* ) WI < llp|lc-ll/ ( ' l lli S .(i + of 


gives the desired upper bound 


\RiN,p,9)\ < Const.\\p\\cA\f {s) \\L™ -p 


P 


(2.17) 


(2.18) 


Here and below Const represents (possibly different) generic constants; also, || ■ \\l™ c indicates the 
L°° norm to be taken over the local support of if p ,o- Note that < oo as long as 9 is chosen 

so that /(•) is free of discontinuities in (x — 9tt,x + 9ir). 


To upperbound the truncation error we use Young’s inequality followed by (2.4), 


\T(N,p,0)\ < ||(iSjv/ — /) * i'tfpfi ~ S n iP p> o)\\l° 


< ||SW/ - /IIl 1 • W’fpfi - SNifp,e\\L°° < M\\SNf - /Hl 1 • \\'f’ P ,e\\c 3 j^z y• (2.19) 


Leibnitz rule yields, 


W$l < r (.+D£ ( * ) | p (.-*)| . |D»| 

k =0 ' ' 


< \\p\\c° 


i +p 


S+l 


( 2 . 20 ) 


and together with (2.19) we arrive at the upper bound 


\T[N ,p, 9; f)\ < Const\\S N f-f \\ L i • ||p|| c . 


(1 +p)N f l+p 


9 


N9 


( 2 . 21 ) 


A slightly tighter estimate is obtained by replacing the L 1 — L°° bounds with L 2 bounds for /’s 
with bounded variation, 


I T(N,p, 9) I < \\S N f - f \\ L 2 X II S N ipp,e - V’p.el \l 2 < 

< Const.\\f\\ B v • N-V* x y$\\ L2 ■ TV-^-V 2 ), 


and (2.20) then yields 


\T(N,p,9)\ < Const.11pile's • N 


(l±p\ 
\N9 J 


s+l 


( 2 . 22 ) 


(2.23) 
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Using this together with (2.IS), we conclude with an error bound of E{N,p,0\ f(x)), 


\ip P ,e * S N f(x ) - f(x )| < Const\\p\\c° 


N 


1 + 


p\ 


s+l 


NO J 


+ P 


ii/ 1 ' 


s )l 


Lf° (x) 

IOC ' > 


Vs, (2.24) 


where ||/( s )||l°° = sn Pye(x-Oir,x+$ir) |/^| measures the local regularity of /. It should be noted 
that one can use different orders of degrees of smoothness, say an r order of smoothness for the 
truncation and s order of smoothness for the regularization, yielding 


\E(N,p, 9; f(x))\ < Const. 


I C' 


■ N 


1 + 


p\ 


r +1 


NO ) 


+ P c« -P 


P 


ii/ ( 


S )| 


Vr, s. 
(2.25) 


2.3 Fourier interpolant - error analysis for pseudospectral mollifier 

The Fourier interpolant of a 27r-periodic function, /(•), is given by 


lNf{y) ■= ^2 Ike 

\k\<N 


iky 


fk ■= 


l 

2 N 


2N—1 


J2 


,-iky v 


(2.26) 


y=0 


We observe that the moments computed in the spectral projection (|2.1|) are replaced here by the 


corresponding trapezoidal rule evaluated at the equidistant nodes y u = v = 0,1,..., 2 N — 1. 
It should be noted that this approximation by the trapezoidal rule converts the Fourier-Galerkin 
projection to a Pseudo Spectral Fourier collocation (interpolation) representation. It is well known 
that the Fourier Interpolant also enjoys spectral convergence, i.e. 


\lNf(x) - f{x )| < Const\\f\\ C s 


1 


Vs. 


(2.27) 


N s ~ l ’ 

Furthermore, if /(•) is analytic with analyticity strip of width 2p, then S]\rf(x) is characterized by 
an exponential convergence rate |Ta94|| 


I Sjv/(x) — f(x) | < Const v ■ Ne Nv . 


(2.28) 


If, however, /(•) experiences a simple jump discontinuity, then the Fourier Interpolant suffers from 
the reduced convergence rate similar to the Fourier projection. To accelerate the slowed convergence 
rate we again make use of our two-parameter mollifier (|2.9| ). When convolving In/{x) by our two 
parameter mollifier we approximate the convolution by the Trapezoidal summation 


iP P ,e*iNf(x) ~ f(y^p,s( x -y^)- 


2N-1 


(2.29) 


u=0 


We note that the summation in ( 2.29|) bypasses the need to compute the pseudo spectral coefficients 
fk. Thus, in contrast to the spectral mollifiers carried out in the Fourier Space [MM078], we are 
able to work directly in the physical space through using the sampling of /(•) at the equidistant 
points, f(y u ). 

The resulting error of our discrete mollification at the fixed point x is given by 
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2N—1 

E(N,p, 9) := ^2 f(yu)^p,e( x ~ Vv) ~ /( x )' ( 2 - 30 ) 

u =0 

As before, we decompose the error into two components 


7r 

N 


E(N,p,0) = 

= : A(N,p,8) + R(N,p,8), 


2N-1 \ 

■ 22 -Vv)- f* ^p,o I + (/ * 4> P ,e - f) 

u=0 ) 


(2.31) 


where R(N,p,8 ) is the familiar regularization error, and A(N,p,6) is the so-called aliasing error 
committed by approximating the convolution integral by a Trapezoidal sum. It can be shown that, 
for any m > 1/2, the aliasing error does not exceed the truncation error, e.g., | Ta94| , (2.2.16)], 


\\A(N,p,0)\\l°° <M m \\T(N,p,e-f^)\\ L oo . N W 2 ~ m \ m> 1/2. 


(2.32) 


We choose m = 1: inserting this into ( |2.21| ) with / replaced by f' , and noting that ||— / IIli < 
Const\\f\\BvVN , we recover the same truncation error bound we had in ( 2.21 ), 


i / -| j \ S+l 

\\A(N,p, 0)\\ L °° <Const.\T(N,p,d-f')\-^=<Const.\\p\\cs-N 2 ^-^j . (2.33) 

Consequently, the error after discrete mollification of the Fourier Interpolant satisfies the same 
bound as the mollified Fourier projection 


\E(N,p,8-J(x))\ < ConstHpHc* 




Vs > 1/2. (2.34) 


We close by noting that the spectral and pseudospectral error bounds, (|2.24| ) and (|2.34j) , are of 
the exact same order. And as before, one can use different orders of degrees of smoothness for the 
regularization and aliasing errors. 


2.4 On the choice of the (9,p) parameters — spectral accuracy. 

We now turn to asses the role of the parameters, 8 and p, based on the spectral and pseudospectral 
error bounds ( 2.24| ) an d (|2.34| ). We first address the localization parameter 8. According to the first 
term on the right of ( |2.24j ), and respectively - ( 2.34 ), the truncation, and respectively - aliasing 
error bounds decrease for increasing 0’s. Thus we are motivated to choose 8 as large as possible. 
However, the silent dependence on 8 of the regularization error term in ( [2.24 ) and Q2.34| ) appears 
through the requirement of localized regularity, i.e. = sup ye ^ x _ dwx+gn )\f^ s \y)\ < oo. 

Hence, if d(x) denotes the distance from x to the nearest jump discontinuity of /, 


d(x) := dist(x, sing supp /), 


(2.35) 


we then set 


8 := 



7r 


(2.36) 
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This choice of 9 provides us with the largest admissible support of the modifier ip Pi g, so that 
i/jp t g* f(x) incorporates only the (largest) smooth neighborhood around x. This results in an adap¬ 
tive modifier which amounts to a symmetric windowed filter of maximal width, 2 d(x), to be carried 
out in the physical space. We highlight the fact that this choice of an x-dependent, 6{x) = d(x)/Tr, 
results in a spectral modifier that is not translation invariant. Consequently, utilizing such an 
adaptive modifier is quite natural in the physical space, and although possible, it is not well suited 
for convolution in the frequency space. 

How can we find the nearest discontinuity? we refer the reader to [GeTa99, GeTaOOa] , |GeTa00b| , 
where a general procedure to detect the edges in piecewise smooth data from its (pseudo-)spectral 
content. The procedure - carried in the physical space, is based on appropriate choice of concentra¬ 
tion factors which lead to (generalized) conjugate sums which tend to concentrate in the vicinity of 
edges and are vanishing elsewhere. The locations (and amplitudes) of all the discontinuous jumps 
are found in one global sweep. Equipped with these locations, it is a straightforward matter to 
evaluate, at every x, the appropriate spectral parameter, 6{x) = d(x)/ir. 

Next we address the all important choice of p which controls how closely 'tjj p p possesses near vanish¬ 


ing moments of increasing order, (2.11). Before determining an optimal choice of p let us revisit the 
original approach taken by Gottlieb and Tadmor GoTa85) l . To this end, we first fix an arbitrary 
degree of smoothness s, and focus our attention on the optimal dependence of p solely on N. With 
this in mind, the dominant terms of the error bounds (2.24) and ([2.34 ), are of order ( p/N) s and 
p~ s , respectively. Equilibrating these competing terms gives p = y/N, which results in the spectral 
convergence rate sought for in [ GoTa85|1 , namely, for an arbitrary s 


I E(N,p,t 


0 =Vn — 


< Const s qN s / 2 . 


(2.37) 


Although this estimate yields the desired spectral convergence rate sought for in GoTa85| , it 
suffers as an over-pessimistic restriction since the possible dependence of p on s and 9 were not 
fully exploited. In fact, while the above approach of equilibration with p-depending solely on N 
yields p = N 0,5 , numerical experiments reported back in the original | GoTa85| 1 have shown that 
when treating p as a fixed power of IV, p = N@, superior results are obtained for 0.7 < /3 < 0.9. 
Indeed, the numerical experiments reported in §3.2 below and which motivate the present study, 


clearly indicate that the contributions of the truncation and regularization terms are equilibrated 
when p ~ N. Moreover, the truncation and aliasing error contributions to the error bounds (2.24) 
and ( 2.34|) predict convergence only for x’s which are bounded away from the jump discontinuities 
of /, where 9(x) > p/N. Consequently, with 9(x) := d(x)/-n and /(•) having a discontinuity, say at 
xq, convergence can not be guaranteed in the region 


/ P , P \ 

(X° - — 7T,X 0 + Jf™)- 


(2.38) 


Thus, a non-adaptive choice of p - chosen as a fixed fractional power of N independent of 9(x), 
say p ~ Vn, can lead to a loss of convergence in a large zones of size 0(N 1 / 2 ), around the 
discontinuity. The loss of convergence was confirmed in the numerical experiments reported in 
S |3.2| . This should be contrasted with the adaptive mollifiers introduced in the next §|3], which 
will enable us to achieve exponential accuracy up to the immediate, 0(1/N) vicinity of these 
discontinuities. We now turn to determine an optimal choice of p by incorporating both - the 


distance to the nearest discontinuity, d(x), and by exploiting the fact that the error bounds (|2.24[) 


and (2.34) allow us to use a variable degree of smoothness, s. 
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3 Adaptive Modifiers — Exponential Accuracy 


Epilogue - Gevrey regularity. The spectral decay estimates (|2.2| ) and ( 2. 27] ) tell us that for 
Cq° data, the (pseudo) spectral errors decay faster than any fixed polynomial order. To quantify 
the actual error decay, we need to classify the specific order of Cq° regularity. The Gevrey class, 
G a ,a > 1, consists of p’s with constants p := rj p and M := M p , such that the following estimate 
holds, 


sup|p (s) (x)| < s = 1,2,... (3.1) 

ieS V 

We have two prototypical examples in mind. 

Example 1 A bounded analytic function p belongs to Gi with M p = sup^gjj |p(x)| and 2p p equals 
the width of p’s analyticity strip. 

Example 2 Consider a Cq°(— 7r, it) cut-off function depending on an arbitrary constant c > 0, 
which takes the form 


Pc{x) = < 




0 


x\ > 7T 


(3.2) 


In this particular case there exists a constant A = A c such that the higher derivatives are 
upper bounded by] 3 4 ! 


/ cx 2 \ 

s = 1,2,... (3.3) 

The maximal value of the upper bound on the right hand side of (|h^) is obtained at x = x m ax 
where x 2 nax — ir 2 ~ —7r 2 c/s; consult]]. This implies that our cut-off function p c admits G 2 
regularity, namely, there exists a constant rj c := \ c ir 2 c such that 

sup |pP s V.t)| < Const c ■ s! ( —^ e~ s < Const c ^-^— s = 1,2,... (3.4) 

rrSSR \VcJ Vc 

We now turn to examine the actual decay rate of Fourier projections, 15/vp — p|, for arbitrary 
G Q -functions. According to (^) combined with the growth of ||p||c s dictated by (|3.1|), the L°° 
error in spectral projection of a G a function, p, is governed by 

\Snp(x) - p(x)] < Const. N ^ , s = l,2,... (3.5) 

The expression of the type encountered on the right of (|3.5|), (sl) a (r]N)~ s , attains its minimum at 

Smin = (plV) 1/Q , 

3 To this end note that p c (x) = e+(x)e-(x ) with e±{x) := exp(cx/(x ± w) for x £ (—n,n). The functions e±(x) 
upper bounded by |e(^(a;)| < M±s\(\ c \x ± 7r|) _s e±(a:), with appropriate A = A c , p. 73] 

4 For large values of s, the function 1 0 .( 33 )| —s • exp(aa(x) + (3/a(x)) with fixed a and /3 is maximized at x = x m ax 
such that a(xmax) ~ —/3/s. I 11 our case, a(x) — x 2 — n 2 and /3 ~ cn 2 
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■ («!)“ 

mm 7 ——— ~ mm . 

s ( r]N) s s yr)e a N 


ex. \ s 

,s ' g— a (77AT ) 1 / 0 


(3.6) 


Thus, minimizing the upper-bound in 
of fractional order 


at s — Sr 


= (r/A^) 1 /", yields the exponential accuracy 


\Snp(x) — p{x)\ < Const ■ Ne a ( r i N ) 1/a ^ p £ G 0 


(3.7) 


The case a = 1 recovers the exponential decay for analytic p’s, (|2.5|), whereas for a > 1 we have 


exponential decay of fractional order. For example, our G 2 cut-off function p = p c in (|3.2[) satisfies 
m with ( 77 , a) = (rf c , 2), yielding 


\SpPc(x) — p c (x) | < Const.p ■ e (3.8) 

Equipped with these estimates we now revisit the error decay of spectral mollifiers based on G a 
cut-off functions p. Both contributions to the error in ( 2Tj) — the regularization R(N,p,0), and 
the truncation T(N,p,9 ) (as well as aliasing A(N,p,9 ) in ( j2.31 )), are controlled by the decay rate 
of Fourier projections. 


3.1 The (6,p) parameters revisited — exponential accuracy. 

We assume that /(•) is piecewise analytic. For each fixed x, our choice of 8 = Q(x) = d(x)/ir 
guarantees that f(x — 8y) is analytic in the range |y| < 7 r and hence its product with the G a (— tt, 7 r) 
function p{y) yields the G a regularity of g x (y) = f(x — 9y)p{y) — f(x). According to s , the 
regularization error, R(N,p,9 ) is controlled by the Fourier projection of g x (-), and in view of its 


G a regularity, (|3.7|) yields 


\R(N,p,9)\ = \(S p g x (y) - g x (y))\ y=0 \ <Const p -pe a (vp) 1/a 
For example, if p = p c we get 

| R Pc (N,p, 0)| < Const c ■ p ■ 


(3.9) 


(3.10) 


Remark. It is here that we use the normalization, p(0) — 1 = g x (y = 0) = 0, and (3.£) shows that 
one can slightly relax this normalization within the specified error bound 


|p(0) - 1| < Const. e~ Q ^ 1/o 


(3.11) 


Next we turn to the truncation error, T(N,p,9). According to (|2.19[) , its decay is controlled by 
the Fourier projection of the localized Dirichlet kernel 'ipp.o(x) = (f)- Here we shall need the 

specific structure of the localizer p(x) = p c (x) in (5]^). Leibnitz rule and (^) yield 


ivPmi < E 


k=0 


k , U k \x)\ ■ \D<f~ k Xx)\ 


< Const, s! I ^2 


P S k / I 2 21 \ —fc | \x 2 -7T 2 

(Vc\X -TT I) • e v 


\k =0 


(s-k) 


cx 


< Const. 


SI 


/ \ I 9 01 ' 


p\ c \x 2 — 7T 2 | + 


cx 


X 2 — 7T 2 




















12 


E. Tadmor and J. Tanner 


which after dilation satisfies 




• e. V 


( p\ c \a(x)\ cx 2 

9 2 a(x ) 


a(x) := x l — tx 2 9 2 {x) 


(3.12) 


Following a similar manipulation we used earlier, the upper bound on right hand side of the (3.12) 
is maximized at x = x ma x with x^ ax — vr 2 0 2 ~ —ck 2 9 2 /s 1 which leads to the G 2 -regularity bound 
for (where as before, p c := A c 7r 2 c) 


sup \ipp S l{x)\ < Const ■ s! 
a;S5R ’ 


VcOe 


e P11c/s < C osn l' 


( s Q 2 

(VcP) S 


s PVc/s 


S = 1 , 2 ,... 


With (3.13) we utilize (2.22) to obtain the following precise bound of the truncation error 

\T{N^6)\<Const.j^y s e^l s 


(3.13) 


(3.14) 


To minimize the upperbound ( 3.14 ), we first seek the minimizer for the order of smoothness, s = s p , 
and then optimize the free spectral parameter p < N for both the truncation and regularizations 
errors. We begin by noting that a general expression of the type encountered on the right of (|3.14), 


(-0 


n 2 


0 PVc/s 


(■ ricON) s 

is minimized at the p-dependent index s r 




such that 


d s (log M{s,p)) | s=Smin = log( 


pile 


VcON' 


= 0 


(3.15) 


Though we cannot find its explicit solution, (3.15) yields a rather precise bound on s m in which turns 
out to be essentially independent of p. Indeed, for the first expression on the right of ( 3.15| ) to be 
positive we need s m i n = \fJJrp£N with some f3 > 1. Plugging this expression, s m i n = \/p7j c 9N, 
into ( 3. 15| ) , we find that for p < N we must have, log/3 = log(sf nin /r] c 9N) = p c p/s^ lin < 1/(59. We 
therefore set s ~ s m i n of the form 


s = sj (5r] c 9N, 1<(5< 1.764, 

so that the free (3 parameter satisfies the above constraint^ /31og/3 < 1 < 1/9. The corresponding 
optimal parameter p is then given by 


Pmin = — ■ (log -^rr)| s=s min = « ' 0 < K =: (3 log (3 < 1(< ^). (3.16) 

Pc PcP^ V 

We conclude with an optimal choice of p of order 0(9N), replacing the previous choice, ( P7D , of 
order 0(\^N). The resulting exponentially small truncation error bound, (3.14), now reads 


\T(N,p,9)\ < Const.j^^e ^, s= ^ 




/ p\2V0Vc8N 

w 


1 < (3 < 1.764. 


(3.17) 


'(PcONY 

With this choice of p = p m i n in ( |3.16 ) we find essentially the same exponentially small bound on 
the regularization error in (3.10), 


5 R.ecall that 9 = 9(x) := d(x)/n < 1. 
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/ i \ 2 V/3 log 0 - r]cON 

\R(N,p,9)\ < Const- ON (-J . (3.18) 

Figures EOI(f) and [h^(f) below confirm that the contributions of the truncation and regularization 
parts of the error are of the same exponentially small order up to the vicinity of the discontinuous 
jumps with this choice of optimal p ~ Nd(x)/ir , in contrast to previous choices ofp = 0(/V 7 ), 7 < 1, 
consult Figures |3d|(b)-(d) and ^|(b)-(d). 

We summarize what we have shown in the following theorem. 


Theorem 3.1 Given the Fourier projection, <S)v/(-) of a piecewise analytic /(•), we consider the 
2-parameter family of spectral mollifiers 


■= ^c(|)A>(f )> 


Pc ■= e 


cx 


x 2 — vr 2 


c > 0 , 


and we set 


= 9{x) := 


d{x) 


7T 


d(x) = distfx , sing suppf) 


p = p(x) ~ k ■ 6{x)N, 0 < k = /3 log < 1. 


(3.19) 

(3.20) 


Then there exist constants, Const c and p c , depending solely on the analytic behavior of /(•) in 
the neighborhood of x, such that we can recover the intermediate values of f(x) with the following 
exponential accuracy 


\ip P ,0 * S N f(x) - f(x) | < Const c ■ ON 


p \ 2 yJKrj c 0 ( x)N 


l</3< 1.764. 


(3.21) 


Remark. Theorem |3.1| indicates an optimal choice for the spectral modifier, 'fp.o, based on an 
adaptive degree of order p = k6(x)N, with an arbitrary free parameter, 0 < n = /? log /3 < 1. We 
could further optimize the error bound ( |3 . 21 ) over all possible choices of (3, by equilibrating the 
leading term in the truncation and regularization error bounds so that 

p\ 2 VP-VoO{x)N / 1 \ 2 Vl 3Io S P-Vc9(x)N 


— 


with the minimal value found at log/l* = (3 — \/5)/2, the corresponding k* := /3* log/3* = 0.5596, 
and k^/tt = 0.8445 leading to an error bound, 


| il>p,e * S N f(x ) - f(x) | < Const c ■ d(x)N 


0.845-y/ r / c d ( x)N 


(3.22) 


Although the last estimate serves only as an upperbound for the error, it is still remarkable that the 
(close to) optimal parameterization of the adaptive modifier is found to be essentially independent 
of the properties of /(•). 

Similar result holds in the pseudospectral case. In this case, we are given the Fourier interpolant, 
In f{x) and the corresponding discrete convolution is carried out in the physical space with overall 
error, E(N,p,6; f) = ifpfi * In f(x) — f(x), which consists of aliasing and regularization errors, 
( 2.31 ). According to ( 2.32| ) , the former is upper bounded by the truncation of f , which retains the 
same analyticity properties as / does. We conclude 


















14 


E. Tadmor and J. Tanner 


Theorem 3.2 Given the equidistant gridvalues, {f (x v )}o<„< 2 N-i of a piecewise analytic /(■), we 
consider the 2-parameter family of spectral mollifiers 


and we set 


4>p,o( x ) := ]jPc(]j) D p (§)> 


O' 


Pc := e 


cx 


c > 0, 


d(x\ 

6 = 0(x) := -, d(x) = dist(x,sing suppf) 

7r 

p = p(x) ~ k • 6(x)N, 0 < k = piog/3 < 1. 


(3.23) 

(3.24) 


Then, there exist constants, Const c and p c , depending solely on the analytic behavior of /(•) in 
the neighborhood of x, such that we can recover the intermediate values of f(x) with the following 
exponential accuracy 


7r 

N 


2N—1 


^2 4’ P ,e(x - y v )f{y v ) - f(x) < Const c ■ (d{x)Nf 


u=0 


p\ 2y/0 Vc9(x)N 


1 < p < 1.764. (3.25) 


3.2 Numerical experiments 


The first set of numerical experiments compares our 
involving the same choice of /(•) = /i(-) 


results with those of Gottlieb-Tadmor, [GoTa85||, 


J sin(x'/2) x 
\ — sin(x/2) x 


A second set of results is demonstrated with 



(3.26) 


/2(*) 


(2e 2x - 1 - e 1T )/(e 7r 
— sin(2x/3 — 7 t/3) 



(3.27) 


This is a considerable challenging test problem: fi has a jump discontinuity at x = 7r/2 and due to 
the periodic extention of the Fourier series two more discontinuities at the boundaries x = 0, 27r. 
Moreover, a relatively large gradient is formed for x ~ 7r/2—, and the sharp peak on the left of 
x = 7t/ 2 is met by a jump discontinuity on the right. 

For the computations below we utilize the same localizer p = p c as in ([h^), with c = 10. In the 
first case, fi has a simple discontinuity at x = it so the 6 parameter was chosen according to (2.36), 
6 = 9{x) = min(|x|, \x — 7r|)/7r. In the second case of f 2 {x) we set 


6{x) = [min(x, 7r/2 — x)+ + min(x — 7 t/2, 2ir — x)+] /it. 
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Since the error deteriorates in the immediate vicinity of the discontinuities where 9{x)N ~ 1, a 
window of minimum width of 9 m i n = min{0(x), 1/4 A} was imposed around Xq = {0, 7t/2, 7t, 27t}. 
More about the treatment in the immediate vicinity of the discontinuity is found in §4.1. 

The different policies for choosing the parameter p are outlined below. In particular, for the near 
optimal choice recommended in Theorems 3.1 and 3.2 we use a modifier of an adaptive degree 
p = k6(x)N with k = 1 /yfe = 0.6095 ~ re*. 

We begin with the results based on the spectral projections, S'jv/i and Sjv/2- For comparison 
purposes, the exact convolution integral, 1 / 7.0 *5/*// was computed with composite Simpson method 
using points, and the mollified results are recorded at the left half points, v = 0 , 

1,..., 149. Figure [O] shows the result of treating the spectral projection 'tp p p * Sjv/i based on 
A = 128 modes, for different choices of p’s. Figures [h^ show the same results for f 2 (x). It is 
evident from these figures, figure |3T] (e)-(f) and figure |3.2| (e)-(f), that best results are obtained with 
p = 9(x)N/y/e, in agreement with our analysis for the optimal choice of exponentially accurate 
modifier in (|3.22|) . We note that other choices for p ~ A 7 , lead to large intervals where exponential 
accuracy is lost due to the imbalance between the truncation and regularization errors, consult 
cases (a)-(d) in figures |3T[ - |3.2| . As we noted earlier in ( 2.38 ), a nonadaptive choice of p independent 
of 9(x) leads to deterioration of the accuracy in an increasing region of size ~ pn/N around the 
discontinuity, and the predicted locations of these values, given in table 3J, could be observed in 
figures [ 0 (a)- (d) for the function fi(x). 


p\N 

32 

64 

128 

C 

oc 

1.6 

1.8 

2.0 

A u - 5 

2.6 

2.7 

2.9 

A a2 

2.9 

3.0 

3.1 


Table 3.1: Predicted location where spectral convergence is lost at \x — xq\ ~ pir/N. 


Figure [h^ illustrates the spectral convergence as A doubles from 32 to 64, then to 128H. The 
exponential convergence of the near optimal adaptive p = 9(x)N/yfe can be seen in figure £F3|(e)- 
(f), where the log-slopes are constants with respect to d(x ) (for fixed A) and with respect to A 
(for fixed x). 

Next, the numerical experiments are repeated for the discrete case, using discrete mollification of 
the Fourier interpolant. Given the gridvalues of fi(x u ) and ^(xj,) at the equidistant gridpoints 
x v = un/N, we recover the pointvalues at the intermediate gridpoints f(x u . 1 - 1 / 2 )• A minimal window 
width of 9 m i n = min(0(x), 2n/N) was imposed in the immediate vicinity of the discontinuities to 
maintain a minimum number of two sampling points to be used in the discrete mollification. 
Compared with the previous mollified results of the spectral projections, there are two noticeable 
changes, both involving the non-optimal choice of p ~ A 7 with 7 < 1: (i) The location where 
spectral/exponential convergence is lost is noticeably closer to the discontinuities compared with 
the mollified spectral projections, but at the same time (ii) Much larger oscillations are observed 
in the regions where spectral convergence is lost. Comparing figures [hl| vs. 3.4, and 3.2 vs. 3.5, 
gives a visual comparison for both changes from spectral to pseudo spectral. The deterioration for 
p = A 0,8 and A 0,5 are very noticeable. 


Machine truncation error is at ~16. 






















Figure 3.1: Recovery of fi(x) from its first N = 128 Fourier modes, on 
responding regularization errors (dashed) and truncation errors (solid) o 
spectral modifier based on various choices of p: (a)-(b) p = N 0,5 , (c 
p = Nd(x)/iry/e. 





Figure 3.2: Recovery of /^(x) from its first N = 128 Fourier modes, on 
responding regularization errors (dashed) and truncation errors (solid) o 
spectral modifier based on various choices of p: (a)-(b) p = N 0,5 , (c 
p = Nd{x)/Ti\fe. 
















Figure 3.4: Recovery of 
responding regularizatio: 
spectral modifier ip Pj Q be 
p = Nd(x)/ir\/e. 













Figure 3.5: Recovery of 

j 2 \x) from its N = 

128 eqi 

responding regularizatio 

n errors (dashed) an 

d trum 

spectral modifier ip p p b; 
p = Nd(x)/ir\/e. 

rsed on various choic 

es of p: 
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4 Adaptive Mollifiers — Normalization 


The essence of the 2-parameter spectral mollifier discussed in §|3[ if p fi(x), is adaptivity: it is based 
on a Dirichlet kernel of a variable degree, p ~ 0(x)N, which is adapted by taking into account the 
location of x relative to its nearest singularity, 9{x) = d(x)/ir. The resulting error estimate tells us 
that there exist constants, Const , 7 and a > 1 such that one can recover a piecewise analytic f(x) 
from its spectral or pseudospectral projections, P/v/( - ), 

I tp P ,6*PNf(x) ~ f(x) | < Const. e -(.i d (*)) N ) 1/a ' 

The error bound on the right shows that the adaptive mollifier is exponentially accurate for all x’s, 
except for what we refer to as the immediate vicinity of the jump discontinuities of /, namely, those 
x's where d(x) ~ 1/A. This should be compared with previous, non-adaptive choices for choosing 
the degree of for example, with p ~ y/N we found a loss of exponential accuracy in a zone 
of size ~ 1/y/N around the discontinuities of /. Put differently, there are 0(VN) ’cells’ which are 
not accurately recovered in this case. In contrast, our adaptive mollifier is exponentially accurate 
at all but finitely many cells near the jump discontinuities. According to the error estimates ( |3 .21 ), 
(3.22) and ( |3.25| ), convergence may fail in these cells inside the immediate vicinity of sing suppf, 
and indeed, spurious oscillations could be noticed in figures [hi] and 3.2 . In this section we address 
the question of convergence up to the jump discontinuities. 

One possible approach is to retain a uniform exponential accuracy up to the jump discontinuities. 


Such an approach, developed by Gottlieb, Shu, Gelb and their co-workers is surveyed in ||GoSh95 , 
GoSh98|1 . It is based on Gegenbauer expansions of degree A ~ N. Exponential accuracy is retained 
uniformly throughout each interval of smoothness of the piecewise analytic /. The computational 
of the high order Gegenbauer coefficients, however, is numerically sensitive and the parameters 
involved need to be properly tuned in order to avoid triggering of instabilities, |]Ge97 , GcOO]. 

Here we proceed with another approach where we retain a variable order of accuracy near the jump 
discontinuities, of order 0((d(a:)) r+1 ). Comparing this polynomial error bound against the interior 
exponential error bounds, say 


Const.d{x)N ■ e -o.845 y/ nc d(x)N > (d(x)) r+1 , 


we find that there are only finitely many cells in which the error 
two, is dominant by polynomial accuracy 


d{x) < Const. 


r 2 (logd(x )) 5 

r) c N 


N 


dictated by the smaller of the 


(4.1) 


In this approach, the variable order of accuracy suggested by (4.1), r ~ yj d(x)N , is increasing 
together with the increasing distance away from the jumps, or more precisely - together with the 
number of cells away from the discontinuities, which is consistent with the adaptive nature of our 
exponentially accurate mollifier away from the immediate vicinity of these jumps. The current 
approach of variable order of accuracy which adapted to the distance from the jump discontinu¬ 
ities, is reminiscent of the Essentially Non Oscillatory (ENO) piecewise polynomial reconstruction 


employed in the context of nonlinear conservation laws [ HEOC85 , Sh97). 

How to enforce that our adaptive mollifiers are polynomial accurate in the immediate vicinity of 
jump discontinuities? as we argued earlier in (| 2 . 11 ), the adaptive mollifier fi p y admits spectrally 


small moments of order p s ~ ( d(x)N ) s . More precisely, using (3.7) we find for p € G, 
































22 


E. Tadmor and J. Tanner 


= <5 s q + Const.9 s ■ e p~d(x)N. 


Consequently, C P ,0 possesses exponentially small moments at all x's except for the immediate 
vicinity of the jumps where p ~ d(x)N ~ 1, the same 0(1/N) neighborhoods where the previous 
exponential error bounds fail. This is illustrated in the numerical experiments exhibited in §3.2 
which show the blurring in symmetric intervals with width ~ 1/N around each discontinuity. To 
remove this blurring, we will impose a novel normalization so that finitely many moments of (the 
projection of) ip p .o precisely vanish. As we shall see below, this will regain a polynomial convergence 
rate of the corresponding finite order r. We have seen that the general adaptivity (44) requires 
y 7 d(x)N ; in practice, enforcing a fixed number of vanishing moments, r ~ 2,3 will suffice. 


f rsj 


4.1 Spectral normalization - adaptive modifiers in the vicinity of jumps 

Rather than ipp,e possessing a fixed number of vanishing moments as in standard mollification ( p.7| ), 
we require that its spectral projection, S^ippp, posses a unit mass and, say r vanishing moments, 



y s (S N 'ip p p)(y)dy = 6 s0 


s = 0,1,... ,r. 


(4.2) 


It then follows that adaptive mollification of the Fourier projection, tppp* Sjyf , recovers the point- 
values of / with the desired polynomial order 0(d(x)) r . Indeed, noting that for each x, the function 
f(x — y) remains smooth in the neighborhood |y| < nO = d{x) we find, utilizing the symmetry of 
the spectral projection, / ( S N f)g = f f{S N g), 


/ 7TO 

[f(x ~y)~ f{x)]{S N 'ip p p)(y)dy 

-7 t6 


s =1 


(r + 1)! 
1 \H-1 


-r-/ (s) ( x ) / y s (S N ip> p p)(y)dy + 

. .7—7r 


/ (r+1) (0 / y r+1 {S N ip Pt e){y)dy 


(S N y r+1 )ipp,o(.y)dy < Const. ld(x) + 


N 


The last step follows from an upperbound for the spectral projection of monomials outlined at the 
end of this subsection. 

To enforce the vanishing moments condition ([h^) on the adaptive mollifier, = (p(’)^p('/G))/9, 

we take advantage of the freedom we have in choosing the localizer />(•). We begin by normalizing 


^ P ,e{y) 


•0 P ,e(y) 

fI n *l>P,o(z)dz 


so that i[) Pt Q has a unit mass, and hence (4.2) holds for r = 0, for f SN(i/j Pt g)(y)dy = f tj) p ,g{y)dy = 1. 
We note that the resulting mollifier takes the same form as before, namely 


V’p.e(y) : = l(pcD P )(|), 

where the only difference is associated with the modified localizer, 

1 

Jln‘ l l’p,o(z)dz' 


(4.3) 


Pc(y) = qo ■ p c (y), 


qo = 


(4.4) 
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Observe that in fact, 1/go = f ^ P ,e(z)dz = f ip p (z)dz = (D p *p c )( 0), and that with our choice of 
p = k ■ 9(x)N , we have in view of ([hi]), 


/5c (0) = g 0 = 


(D p * p c ) ( 0) 


= 1 + 0(e), £ ~ d(x)N ■ e 2 C 7 dtf> i p = n-9(x)N 


which is admissible within the same exponentially small error bound we had before, consult ( 3.11 ). 
In other words, we are able to modify the localizer p c (-) —> p c (-) to satisfy the first-order normal¬ 
ization, (|4.2|) with r = 0, while the corresponding mollifier, ( p x D p )q —» ( p c D p )g , retains the same 
overall exponential accuracy. Moreover, using even p’s implies that ip(-) is an even function and 
hence its odd moments vanish. Consequently, ( |4.2| ) holds with r = 1, and we end up with the 
following quadratic error bound in the vicinity of x (compared with (3.22|)) 


\Mx)*S N f(x) - f(x)\ < Const. (d(x) + ^) 2 • e -0.845V^d(*)JV_ 


In a similar manner, we can enforce higher vanishing moments by proper normalization of the 
localizer p(-). There is clearly more than one way to proceed - here is one possibility. In order to 
satisfy (ED with r = 2 we use a pre-factor of the form p c (x) ~ (1 + q 2 X 2 )p c (x). Imposing a unit 
mass and vanishing second moment we may take 

P ,e(y ) = q(pcD p )(|), pc(y) ~ (l + q2y 2 )pc(y), 
with the normalized localizer, p c (y ), given by 

1 + q2y 2 


Pc(y) = 


/!*(! +92(3 Wp,0(2)^ 


Pc(y), 


92 = 


- Jl 7r (^NZ 2 )'i/.’ P) g(z)dz 
j\(SNZ 2 )(^) 2 il) P) e(z)dz' 


(4.5) 


As before, the resulting mollifier i/j Pj q is admissible in the sense of satisfying the normalization ( 3.11 ) 
within the exponentially small error term. Indeed, since f y 2 ^ P ,e(y)dy = (D p -k(y 2 p c (y)))( 0) = Q(e) 
we find 


Pc(0) = 


1 


/!»(! +®(K) 2 )^(y)dy 


= 1 + Const ■ ( d(x)N ) 3 • e 2 V 7 Cp^ p = k ■ 9(x)N. 


1 + 92• s /^ 2 

A straightforward computation shows that the unit mass ippfi has a second vanishing moment 


i: 


y 2 (S N *p p ,g)(y)dy = 


( S N y 2 ) (ao + a 2 (|) 2 ) 'tpp,e(y)dy = 

(S N y 2 )i' P ,e(y)dy + q 2 j (S N y 2 )(^) 2 ip Pt g(y)dy = 0. (4.6) 


Since p c (-) is even, so is the normalized mollifier and hence its third moment vanishes yielding 
a 4th order convergence rate in the immediate vicinity of the jump discontinuities, 

| $ p ,o{x)*S N f)(x) - f(x) I < Const. (d(x) + 1)' 1 • e -o.845Vr MN' 

We close this section with the promised 
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Lemma 4.1 ([Tao]). The following pointwise estimate holds 


\S N {y r )\<{\y\ + ^y. 


To prove this, we use a dyadic decomposition (similar to the Littlewood-Paley construction) to 
split 

= E 2"W2 fc ) 

fc <0 

where ^ is a bump function adapted to the set {tt/4 < \y\ < i r}. 

For 2 k < 1/N, the usual upperbounds of the Dirichlet kernel tell us that 

|5^(^(-/2 fc ))(y)|<2 fc iV/(l + iV|y|). 


Now suppose 2 k > 1/N. In this case we can use the rapid decay of the Fourier transform of i/.’(-/2 k ) 
for frequencies IV to obtain the estimate 


||(l-5 A r)(V’(-/2 fe ))|| 0O <(2 fe IV)- 100 . 


In particular, since supp i/j ~ 1, we have |<Sjv('i/+/2 fe ))(x)| < 1 when |y| 2 k , and \S N (tp(-/2 k ))(y)\ < 

(2 k N)~ 100 otherwise. The desired bound follows by adding together all these estimates over k. 


4.2 Pseudospectral normalization — adaptive modifiers in the vicinity of jumps 

We now turn to the pseudospectral case which will only require evaluations of discrete sums and 
consequently, can be implemented with little increase in computation time. 

Let / * g(x) := f( x ~ yv)g{yv)h denote the (non-commutative) discrete convolution based on 
2 N equidistant gridpoints, y v = vh, h = ir/N. Noting that for each x, the function f(y) remains 
smooth in the neighborhood \x — y\ < tt 8 = d{x), we find 


\Y P ,e * In! {x) - f(x)\ = I ^2^ p ,e(x - y„)[f{y v ) - f(x)\h 


^2 ^p/ (s) (>) ^(x - y v ) s ^ e {x - y u )h + 1 ^ 1 / (r+1) (0 ~ yu) r+1 Y P ,e( x ~ y»)h 


(- 1 ) 


r+l 


s=l 


(r + 1)! 


< Const.(d(x)) 


r+l 


Const 


c, 


r+l 


(r + l)!’ 

provided ip p p has its first r discrete moments vanish, 
2N-1 


^2 ( x ~ yv) s ^p,e{. x - yv) h = S = 0, 1 , 2 , ..., r. 


(4.7) 


v=0 


Observe that unlike the continuous case associated with spectral projections, the discrete constraint 
(4.7) is not translation invariant and hence it requires x-dependent normalizations. The additional 
computational effort is minimal, however, due to the discrete summations which are localized in 
the immediate vicinity of x. Indeed, as a first step we note the validity of (4.7) for x’s which are 
away from the immediate vicinity of the jumps of /. To this end we apply the main exponential 
error estimate ([3. 25) for /(•) = (x — -) s (for arbitrary fixed x), to obtain 
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27V-1 

^2 (x - yu) s i>p,e{x - Vv)h = (x - y)\ y=x + 0(e) = 

u=0 

= 5 s0 + O(£), e~ {d(x)N) 2 . e -\/Const. d {x)- N ' ( 4 8 ) 

Thus, ( [4.7] ) holds modulo exponentially small error for those re’s which are away from the jumps 
of /, where d(x) 7§> 1/N. The issue now is to enforce discrete vanishing moments on the adaptive 
modifier Y P ( X ) = p( x )D p (x) in the vicinity of these jumps, and to this end we take advantage of 
the freedom we have in choosing the localizer p(-). We begin by normalizing 


Y P ,e(y) 


Y P ,e(y) 

Y2=o^Y’pAx ~ y») h 


so that Yp.o( x — •) has a (discrete) unit mass, i.e., ( |4.7|) holds with r = 0. We note that the resulting 
modifier takes the same form as before, namely 


YpAv) ■= \(.PcDp )(|), (4.9) 

and that the only difference is associated with the modified localizer, 


Pc(y) = go ■ p c {y ), go = , — 7 - 77 - ( 4 -i°) 

E^=0 VpA x ~yv) h 

By (|4.8|), the x-dependent normalization factor, qo = qo(x ) is in fact an approximate identity, 


2JV-1 

1/go = E Y P A X - yu)h = 1 + 0(e), £ ~ (d(x)N) 2 • e -Vc° n st.d(x) N , 

u=0 


which shows that the normalized localizer is admissible, |p(0) —1| = |go — 1| < 0 (e), within the same 
exponentially small error bound we had before - consult ( |3.11 ) with our choice of p ~ d(x) ■ N. In 
other words, we are able to modify the localizer p c (-) —* p c (-) to satisfy the first-order normalization, 
(|4.7| ) with r = 0 required near jump discontinuities, while the corresponding modifier, ( p x D p )g —> 
(p c D p )g, retains the same overall exponential accuracy required outside the immediate vicinity of 
these jumps . 

Next, we turn to enforce that first discrete moment vanishes, Ei,( x — yu)Yp,e( x ~ Pv)h = 0, and to 
this end we seek a modified modifier of the form 


YpAy) = 


q{y/8) 


Y.,q^)M x ~y^ 


Y P ,e(y), g(y):=i + giy, 


with q\ is chosen so that the second constraint, ([O]) with r = 1, is satisfied]/] 


gi = 


- yv)Y P A x - yY h 


E: 


(■ x-yv ) 2 

9 


Y P A X -yAh 


(4.11) 


'We note in passing that p c (-) being even implies that ip p ,g(-) is an even function and hence its odd moments 
vanish. It follows that the first discrete moment, ~}2^{x — y^)ip p ,g(x — y v )h vanishes at the gridpoints x = y p , and 
therefore q\ = 0 there. But otherwise, unlike the similar situation with the spectral normalization, q\ ^ 0. The 
discrete summation in q \, however, involves only finitely many neighboring values in the 8 -vicinity of x. 
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Consequently, (4.7) holds with r = 1, and we end up with a quadratic error bound corresponding 

to (p| _ 

$ P ,e * In fix) - fix) | < Const.{d{x)f • e" VConst.d^N _ 


Moreover, (|4.8| ) implies that q\ = O(e) and hence the new normalized localizer is admissible, 
p c { 0) = 1 + 0(e). In a similar manner we can treat higher moments, using normalized localizers, 
Pciy) ~ qiy)pdy) of the form 


^pAv) = \iPc D p)i |), 


' 6 ‘ 


pdy) ■= 


l+qiy + ... q r y r 

Y l v<i(^Tr)' l l ; pA x -yv) h 


pciy)- 


(4.12) 


The r free coefficients of q{y) = 1 + q±y + ...q r y r are chosen so as to enforce (4.7) with the 
first r discrete moments of ip vanish. This leads to a simple r x r Vandermonde system ( - 
outlined in at the end of this section)involving the r gridvalues, {fiy u )}, in the vicinity of x, 
\y v — x\ < 6ix) 7T. With our choice of a symmetric support of size 0(x) = d(x)/7r, there are precisely 
r = 29ir/h = 2IVd(x)/7r such gridpoints in the immediate vicinity of x, which enable us to recover 
the intermediate gridvalues, fix) with an adaptive order (d(x)) r+1 , r ~ Ndix). As before, this 
normalization does not affect the exponential accuracy away from the jump discontinuities, noting 


that p(0) c = 1/qo = 1 + 0(e) in agreement with (|3.1l| ). We summarize by stating 


Theorem 4.1 Given the equidistant gridvalues, {/(x^)}o<i/< 2 Af-i of a piecewise analytic /(•), we 
want to recover the intermediate values fix). To this end, we use the 2-parameter family of pseu¬ 
do spectral mollifiers 

$p,eiy) ■■= ^Pc(|)-Dp(|), P = 0.5596 -ON, c> 0, 
where 6 = 0(x) := d(x)/-7r is the (scaled) distance between x and its nearest jump discontinuity. We 

as the normalizing factor, with 
l + qiy+ ... q r y r 


set pdy) := qiy)e^ v * ^> 1 


qiy) = 


E u qi x -^)M*-y»)h 


so that the first r discrete moments of tbp.oiy) vanish, i.e., 0 holds with r Ndix). 

Then, there exist constants, Const c and t) c , depending solely on the analytic behavior of /(•) in 
the neighborhood of x, such that we can recover the intermediate values of fix) with the following 
exponential accuracy 


2N-1 , s 

jf ^pA x - yv)f{yv) - fi x ) < Constc ■ (d(x )) r+1 y-J 


1 \ 0 . 845 -^/ Vcd{x)N 


r ~ Ndix). (4.13) 


The error bound (4.12) confirms our statement in the introduction of §[|, namely, the adaptivity 
of the spectral modifier in the sense of recovering the gridvalues in the vicinity of the jumps with 
an increasing order, Ndix), proportional to their distance from sing supp f. We have seen that 
the general adaptivity requires r ~ ^Jdix)N, so that in practice, enforcing a fixed number 
of vanishing moments, r ~ 2,3 will suffice as a transition to the exponentially error decay in the 
interior region of smoothness. We highlight the fact that the modified modifier ifpp normalized 
by having finitely many (~ 2,3) vanishing moments can be constructed with little increase in 
computation time and, as we will see in §|4.3| below, it yields greatly improved results near the 
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discontinuities. 

We close this section with a brief outline on the construction of the r-order accurate normalization 
factor q(-). To recover f(x), we seek a r-degree polynomial q(y) := 1 + qiy + ... + q r y r so that (EZD 
holds. We emphasize that the q r ’s depend on the specific point x in the following manner. Setting 
z v := x — y u , then satisfying (|4.7| ) for the higher moments of ip p> g requires 

^2 zl^pAz^h = 0, s = 1,2,... ,r, 

IS 

and with ip Pt g(-) ~ q(-/9)'ip Pt g(-) we end up with 

~ !] ^ P) e( z v)h = - ^2z s v ^ Pi e(z v ), s = 1,2,... ,r. 

IS IS 

Expressed in terms of the discrete moments of i/j, 

a a (z u ) := ^2 (y ) 1+ “V ; p,0(^), a = 1, 2 ... , 2r 

IS 

this amounts to the r x r Vandermonde-like system for {q \,..., q r }, 


aiM 

«2 (z u ) 

* * * Cbr(Zu) 


' qi ’ 



dr+l(Zv) 

dr-\- 2 (^is) 

0^2r {Zis) 


. qr _ 




(4.14) 


Finally we scale q(-) so that ( |4.7| ) holds with s = 0, which led us to the normalized localizer in 

(ED- 


4.3 Numerical Experiments 


Figure [O] (d) the blurring oscillations near the edges when using the non-normalized adaptive 
modifier. To reduce this blurring we will use the normalized if} p g for x’s in the vicinity of the jumps 
where d(x) < 6ir/N. The convolution is computed at the same locations as in section ( |3.2| ), and 
a minimum window width of 9{x) = min(d(a;)/7r, 2-k/N) was imposed, the Trapezoidal rule (with 
spacing of 7r/8000) was used for the numerical integration of {S Ny 2 )^ P fi{y) and ( Sj\ry 2 )y 2 'ip Pt g(y ), 
required for the computation of qo and q 2 in B Figure |4.1| (a)-(d) shows the clear improvement 
near the edges once we utilize the normalized ^ Pt g, while retaining the exponential convergence 
away from these edges is illustrated in figure |4T|(e)-(f). 


We conclude with the pseudospectral case. The 0( 1) error remains in figure 3.4 (d) for the non- 
normalized modifier. The normalization of the discrete modifier in section 4.2 shows that by using 
V’p,0 given in ( |4.12| ), with a 4th degree normalization factor q(-), results in a minimum convergence 
rate of d(.x) 4 in the vicinity of the jumps, and with exponentially increasing order as we move away 
from the jumps. This modification of i() p j) leads to a considerable improvement in the resolution near 
the discontinuity, which could be seen in figure [1 . 2| . Here, normalization was implemented using 
ip Pt 0 in the vicinity of the jumps, for d(x) < 4 tt/N, and the adaptive modifier ip P: g was used for x’s 
’away’ from the jumps d(x) > 47r/V. A minimum window of width 9(x) = min(d(x)/ir, 2ir/N) was 
imposed. 
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Figure 4.1: Recovery of fi(x) (on the left) and /^(x) (on the right) from their N = 128-modes 
spectral projections, using the 4th order normalized modifier (|4.3|) , (|4.5|) of degree p = d(x)N 
Regularization errors (dashed) and truncation errors (solid) are shown on (c)-(d), and Log errors 
based on N = 32, 64, and 128 modes are shown in (e)-(f). 
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(a) 




(b) 



(c) 


(d) 


Figure 4.2: Recovery of fi{x) (a) and foix) (b) from their N = 128-modes spectral projections, 
using the normalized modifier. Log error for recovery of fi(x) (c)and /^(x) (d) from their spectral 
projections based on N = 32, 64, and 128 modes. Here we use the normalized modifier, ifpQ of 
degree p = d{x)N/ 


5 Summary 


In their original work | GoTa85|| , Gottlieb & Tadmor showed how to regain formal spectral con¬ 
vergence in recovering piecewise smooth functions using the 2-parameter family of modifiers 
Our analysis shows that with a proper choice of parameters - in particular, an adaptive choice for 
the degree p ~ d(x)N , hide the overall strength in the method. By incorporating the distance to 
the discontinuities, 9 = d(x)/n along with the optimal value of p, we end up with exponentially 
accurate recovery procedure up to the immediate vicinity of the jump discontinuities. Moreover, 
with a proper local normalization of the spectral modifier, one can further reduce the error in the 
vicinity of these jumps. For the pseudospectral case, the normalization adds little to the overall 
computation time. Overall, this yields a high resolution yet very robust recovery procedure which 
enables one to effectively manipulate pointwise values of piecewise smooth data. 
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